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Abstract: We investigate the measurement of supersymmetric particle masses at the LHC 
in gravitino dark matter (GDM) scenarios where the next-to-lightest supersymmetric part- 
ner (NLSP) is the lighter scalar tau, or stau, and is stable on the scale of a detector. Such 
a massive metastable charged sparticle would have distinctive Time-of-Flight (ToF) and 
energy-loss (dE/dx) signatures. We summarise the documented accuracies expected to be 
achievable with the ATLAS detector in measurements of the stau mass and its momentum 
at the LHC. We then use a fast simulation of an LHC detector to demonstrate techniques 
for reconstructing the cascade decays of supersymmetric particles in GDM scenarios, using 
a parameterisation of the detector response to staus, taus and jets based on full simulation 
results. Supersymmetric pair-production events are selected with high redundancy and 
efficiency, and many valuable measurements can be made starting from stau tracks in the 
detector. We recalibrate the momenta of taus using transverse-momentum balance, and 
use kinematic cuts to select combinations of staus, taus, jets and leptons that exhibit peaks 
in invariant masses that correspond to various heavier sparticle species, with errors often 
comparable with the jet energy scale uncertainty. 



Keywords: 



SUSY, BSM, MSSM 



Contents 



I 



I 



Introduction 

The GDM Benchmark Points 

Fast Monte Carlo Simulation 
|3.l| Parametrisation of Detector 
3.2 Trigger 



s 







3.3| Background 



4.2 



Reconstruction of Masses in GDM Models 

£D] Measuring the Stau Mass 
Recalibrating Taus 
Reconstructing the Decay q 
Reconstructing the Decay q 
Reconstructing v T and %i Decays 



4.3 



4.4 



QXi 
QX2 



grri 
qU- 



4.5 



12 

u 

17 




pSl Conclusions 



22 



1. Introduction 

One of the main motivations for supersymmetry is the existence of a natural candidate 
for cold dark matter (CDM) [1,2] in models in which R-parity is conserved, thus making 
the lightest supersymmetric partner (LSP) stable. Most detector studies of supersym- 
metric models have focused on the constrained minimal supersymmetric standard model 
(CMSSM), in which the soft super symmetry-breaking parameters: the gaugino masses 
mi/2, the scalar masses too and the trilinear parameters Aq, are each assumed to be uni- 
versal at some high scale, and the LSP is assumed to be the lightest neutralino. With the 
precision measurement of the CDM density obtained by combining data from the WMAP 
experiment with other cosmological data [3,4], and taking into account accelerator con- 
straints, the CMSSM parameter space has become quite restricted [5]. Collider signatures 
of these models typically involve missing energy from escaping neutralinos. 

However, there is another plausible CDM candidate in R-conserving supersymmetric 
models, namely the supersymmetric partner of the graviton, the gravitino [6-19]. The grav- 
itino mass, TO3/2, is poorly constrained by accelerator experiments and the astrophysical 
and cosmological constraints are very different from those on the neutralino. In gravitino 
dark matter (GDM) scenarios the next-to-lightest supersymmetric partner (NLSP) has a 
very long lifetime, particularly in gravity-mediated models of supersymmetry breaking, in 
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which the gravitino's couplings to other particles are suppressed by an inverse power of the 
Planck mass. In the case of a neutralino NSLP and an arbitrary gravitino mass, the allowed 
parameter space of CMSSM is expanded and the classic supersymmetric signature of miss- 
ing energy at colliders remains. Here we study instead the intriguing possibility of a charged 
NLSP, in particular the case where the NSLP is the lighter stau, f\. In gravity-mediated 
models with R conservation, the stau would be stable on the scale of a detector, and every 
supersymmetric event would contain a pair of massive metastable charged leptons. Because 
of their large mass and low velocities, these would have distinctive Time-of-Flight (ToF) 
and energy-loss (dE/dx) signatures. Moreover, most supersymmetric events in such stau 
NLSP scenarios also contain a pair of r leptons as well as energetic jets and possibly other 
leptons [20]. 

In Section [2] we review the properties of three benchmark points for GDM models with 
a stau NLSP, which were first proposed in [20] . Subsequently, in Section ||[ we summarise 
the documented expectations for the ATLAS detector response to the metastable f\ with 
limited integrated luminosity, focusing on the accuracy achievable at the LHC in measure- 
ments of the stau mass and momentum for these benchmark points, taking into account 
trigger information and looking at cuts to reject possible backgrounds. In Section |3] we 
use a parametrised fast simulation to investigate the capabilities of ATLAS to measure 
the masses of other supersymmetric particles in GDM models by reconstructing super- 
symmetric cascade decays. For this we use as building blocks the staus themselves, the 
accompanying r leptons, whose full momenta we can reconstruct using transverse momen- 
tum information, hadronic jets and any accompanying charged leptons. In the benchmark 
scenario e, with a relatively light spectrum, many different sparticle species may be re- 
constructed in this way and their masses measured with a high accuracy. In the heavier 
benchmark scenarios £, 77, fewer sparticles can be reconstructed and with worse accuracy, 
even at the fairly high integrated luminosity considered. Finally, in Section |B| we draw our 
main conclusions. 

2. The GDM Benchmark Points 

Three GDM benchmark points with a t\ NLSP, named (e, £, 77), were proposed in [20]. 
These were formulated in the framework of minimal supergravity (mSUGRA) models [21], 
in which the gravitino mass is fixed equal to the universal soft supersymmetry-breaking 
masses of scalar particles at a GUT input scale: 7713/2 = thq, and there is a simple relation 
between the soft trilinear and bilinear parameters: Aq = Bq + rriQ. This relationship can 
be used to fix the ratio of the Higgs vacuum expectation values tan (5 from the electroweak 
vacuum conditions. Further, the value Aq = (3 — v3) tuq found in the Polonyi model 
of supersymmetry breaking in a hidden sector is assumed [22]. The resulting allowed 
GDM parameter space may be displayed as a region in the (777,1/2, mo) plane, subject to 
theoretical, phenomenological and cosmological constraints. 

This region is restricted, in particular, by the requirement of maintaining the cosmolog- 
ical concordance between the values of the baryon-to-entropy ratio inferred from the cosmic 
microwave background radiation and from astrophysical light-element abundances, which 
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constrain NLSP decays and cut the parameter space down to a wedge in the (mi/2,mo) 
plane, throughout which the f\ is the NLSP and is metastable. The shape of this wedge, 
taken from [20], is shown in Fig. |l| as the area restricted to lie below the upper black line. 
The e benchmark point is a collider-friendly point at low values of m\ii and mo sitting 
near the apex of the cosmologically-allowed wedge. It features relatively large production 
cross sections for super symmetric events at the LHC, as seen in Table |l|. The £ and r\ 
points represent more challenging scenarios at higher values of mi/2, and hence with lower 
cross sections. The point £ lies at the boundary of the cosmologically-allowed area, where 
the NLSP lifetime exceeds 10 6 s, whereas r\ is characteristic of models with a relatively 
low t\ lifetime ~ 10 4 s, shown as the lower black line in Fig. |lj, below which additional 
constraints due to hadronic interactions in the early universe become important [25-27]. A 
recent analysis indicates that earlier constraints were overly restrictive, so that the point 77 
lies well within the allowed region [23]. None of these three benchmark points give the full 
amount of CDM required by WMAP, if one takes into account only the gravitino abun- 
dance due to decays of the NLSP after freeze-out. The missing CDM could be provided by 
additional gravitino production mechanisms in the early Universe, or other components of 
CDM such as axions. 
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Table 1: The main supersymmetric production cross sections at the LHC, in pb, including Drell- 
Yan production of f 1 pairs, pair production of strongly- interacting sparticles and associated produc- 
tion of neutralinos and charginos, computed to NLO for the GDM benchmarks using Prospino2 [24]. 
Note that these differ from those given in [20] , mainly because of the larger gluino and squark masses 
obtained from ISAJET (see below) and partly because we use different parton density functions (see 
Section ||) . 

For the purposes of this paper, the effective masses of sparticles at the electroweak 
scale for these benchmark points were calculated by running the universal high-scale masses 
down to low scales using ISAJET 7.69 [28]. This gives somewhat larger masses than found 
in [20] using the SSARD code, in particular for the strongly-interacting sparticles. Both 
the GUT-scale input parameters and the resulting physical masses are given in Table || 
(as in [20], we assume sign(/x) = + and m t = 178 GeV for all three points). The decay 
widths and branching ratios of the supersymmetric particles were re-calculated from these 
masses using the decay code SDECAY 1.1a [29], except for the three-body decays of the 
first- and second-generation right-handed sleptons, which are not included in SDECAY, and 
whose decays instead are taken from ISAJET. 

One of the most important decays of these GDM models, at least for those within the 
reach of the LHC, is that of the right-handed squark: q~R — » qxi with branching ratios of 
almost 100% at all three points. This is followed by the decay x? ~~ > tt\ i n a large fraction 
of events (92/75/69 % for the three points), so that the decays of two right-handed squarks 
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result in the relatively clean final states qqrr f\f\ in a large fraction of the events. From the 
branching ratio of this decay chain shown in the left plot of Fig. [|, we see that this signature 
is a quite generic feature for the class of GDM models with a f NLSP. In the experiments it 
should be possible to trigger on the pair of high-p-p jets produced in the primary qji — > qx\ 
or the pair of tau jets, and one can then identify the stau by its peculiar signature as a 
slow moving muon and measure its mass and momentum to high precision, as discussed in 
the next Section. One can also identify the decay products of the r and reconstruct the 
full r momenta using transverse momentum balance, as we show in Section ||. Thus, these 
events can be fully reconstructed and we can make a precise determination of the masses 
of all three supersymmetric particles involved in the decay cascade. From this decay chain 
one may also be able think of reconstructing the gluino g from its decays into squark-quark 
pairs, as we discuss later for benchmark point e. 

The decays of left-handed squarks have in general more leptons in the final state, 
for example when decaying via X2 ~* an d 0-l — ► ^Xi> with the total branching ratio 
shown in the right plot of Fig. |l[ With knowledge of the Xi mass, this makes possible the 
reconstruction of the left-handed slepton and the x.2- However, left-handed squark decays 
into charginos are more difficult to use because of neutrinos in the final state, either directly 
from the chargino decay, or from the leptonic decay of a W. 

Whilst the branching ratio of the left-handed squark is too low to search for events 
with two such decay chains, the considerations for triggering, particle identification and 
reconstruction are similar to those for the right-handed squark. The appearance of hard 
leptons help identify the decay chain, but give additional combinatorial difficulties in par- 
ticle identification as we will see in Section Nevertheless, it seems possible to reconstruct 
this decay chain, at least for the benchmark point e, as we discuss later. 




Figure 1: Total branching ratios for the qn — > qx\ — > qrf\ (left) and q~L — > qx® —> Q^l — > Q^Xi ~* 
qiirfi (right) decay chains in the (nT,i/2) m o) parameter plane. The solid black lines indicate the 
boundary of the cosmologically allowed region of parameter space explored in [20] . 
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At the e benchmark point the super symmetric events are dominated by gluino/squark 
production followed by cascade decays. For the low cross-section benchmark points £ and 
rj, we see from Table [l] that these events make up a smaller fraction of the total number of 
supersymmetric events. Here Drell-Yan and associated neutralino and chargino production 
become more important. These events feature less jet activity and heavier, slower staus, 
and may lead to some trigger problems (see Section p\2| ) . By considering only the Drell-Yan 
cross-section, taking into consideration triggering and reconstruction efficiencies for stau 
pair production, one could investigate the reach of the LHC in setting model-independent 
exclusion limits on the mass of the stau. However, this would require a very detailed 



understanding of the effects of the background cuts suggested in Section [jyj on SM events 
in the limit where the background tends to zero, and is perhaps best left to a full detector 
simulation. Here we will instead focus on some interesting decay chains starting from the 
associated production. The branching ratio of the decay xf — > tv t — > tWt\ is shown in 
Fig. |2[ and is almost identical to that for the decay \% ~^ v T v T — » v T Wfi, as the decays 
differ essentially only by the mass of the tau. If we can reconstruct the W produced in 
both cases from its hadronic decay and combine with the correct stau candidate, this gives 
us direct access to the tau-sneutrino and, possibly, chargino masses. 

Among the lighter particles of Table [2] that one would expect to be copiously produced 
at the LHC in the case the e benchmark, we have not yet discussed decays involving the 
£r, the right-handed selectron or smuon. These are dominantly produced in decays of 
neutralinos, the only exception being the direct pair-production of the sleptons, which 
have low cross sections. In fact, only in the decay of the lightest neutralino is there a 
significant rate, due to the preference for decays to left-handed sleptons, and even here it is 
suppressed relative to decays to the lighter stau. The £r decays in turn almost exclusively 
to the three bodies £rf\. In the right plot of Fig. we show the total branching ratio of 
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Figure 2: Branching ratio for the decay chains xf — > tv t — > tWt\ (left) and x\ — * ££r — > ££tt\ 
(right) in the (ra^2,mo) parameter plane. The solid black lines indicate the cosmologically allowed 
region of parameter space explored in [20]. 
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the decay Xi ~* ~> Utt~\ . One could hope to reconstruct the £r mass from this decay, 
starting from the stau. However, the low branching ratio at the point e, combined with 
the soft spectra of the leptons from the decay, due to the small mass differences, makes 
this very difficult. Our attempts leave us with only a few events. The larger branching 
ratios and harder leptons of the two other benchmarks are unfortunately balanced out by 
the overall much lower cross sections. It thus seems to be very challenging to reconstruct 
the £r in GDM models at the LHC, and we do not discuss the £r further. 

The first four decay chains we discussed above will be the main focus of our at- 
tempts at reconstructing supersymmetric masses in Section f|. The fact that, as shown 
in Figs. H and [2|, the total branching ratios for these decay chains are relatively high across 
large regions of the (^1/21^0) plane, with a dominating branching ratio for the chain ini- 
tiated by qR — > qx± in all the explored range of the allowed parameter space, up to values 
of mi/2 ~ 2 TeV, suggests that the techniques discussed here would be of wide utility at 
the LHC. 

For the most part, we limit our discussions to simulations of the high-cross-section 
benchmark point e, and only summarise our results for the other two points. We draw 
attention to instances where there are important differences between the benchmark points, 
and we refer the reader to [20] for further details about them. 

3. Fast Monte Carlo Simulation 

For the fast simulation we have generated proton-proton collisions at the LHC energy using 
PYTHIA6.326 [30] and CTEQ 5M1 parton distribution functions [31]. For the e benchmark 
point we show results on mass measurements for a total integrated luminosity of 30 fb _1 , 
corresponding to the planned initial running of the LHC at low luminosity. In the cases of 
the C and r\ points, which have relatively low cross sections, we show results for statistics 
equivalent to 300 fb _1 . These may then be viewed as estimates of the ultimate precisions 
achievable for these benchmark points without an upgrade of the LHC luminosity. 

3.1 Parametrisation of Detector 

The detector simulation has been carried out using the generic LHC detector simulation 
AcerDETl.O [32], taking the ATLAS detector as our model for a LHC detector. Here 
we give a short summary of the most important choices made for the AcerDET settings: 
We consider a lepton to be identified if it has transverse momentum pt > 5(6) GeV 
and I T]\ < 2.5 for electrons (muons). A lepton is considered isolated if it is at a distance 
AR > 0.4, where AR = W (Acp) 2 + (Arj) 2 , from other leptons and jets, and if the transverse 
energy deposited in a cone AR = 0.2 around the lepton is less than 10 GeV. Jets are 
reconstructed by a cone-based algorithm from clusters and are accepted if the jet has 
Pt > 15 GeV within a cone AR = 0.4. The jets are re-calibrated using an included flavour- 
independent parametrisation, optimised to give the correct resonance mass in dijet decays. 
We use the px parametrised 6-tagging efficiency and light jet rejection, for a low-luminosity 
environment, given in [33]. For tau jets we use a lower cut on transverse momentum of 
Pt > 20 GeV, and the parametrised tau-tagging efficiency and rejection factors from [34]. 
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For the detector's response to staus, we take our 
parametrisation of the momentum resolution from [35, 
36], which performed a full simulation of the response 
of the ATLAS muon system to a metastable stau of 
mass 101 GeV, in a Gauge Mediated Supersymmetry 
Breaking (GMSB) model. This momentum dependent 
resolution is used to smear the stau momentum in the 
detector simulation. We have performed a smearing 
with width a p , given by 



= kip®k 2 \ 1 + 



2 

m~ 
n 



h 
p ' 



(3.1) 



where k\ is the parameter of the sagitta measurement 
error, k 2 represents the multiple scattering and k% the 
fluctuation of energy loss in the calorimeter. From [35] 
we take the values k\ = 0.0118, k 2 = 0.02 and k^ = 89. 

For the measurement of the velocity of staus and 
muons, we use the parametrisation of the velocity res- 
olution from [35,36]. In our simulation we smear the 
true velocity by: 



= 0.028/3, 



(3.2) 



which was found by looking at the fit quality of tracks 
in the muon system as a function of the particle's 
assumed arrival time in the Monitored Drift Tubes 
(MDT), subsequently minimising the x 2 as a function 
of the arrival time to find the correct Time-of-Flight 
(ToF). This is a conservative estimate of the resolu- 
tion as the simulation was carried out in the central 
detector, for rj = 0.1. For larger pseudo-rapidities 
the resolution is expected to improve, due to a longer 
flight path. An alternative way to measure the ToF 
would be to use the timing of hits in the Resistive 
Plate Chambers (RPC), responsible for muon trigger- Table 2: Proposed GDM bench- 
ing, which may be able to provide similar precision. mark points taken from [20]. Mass 

spectra [GeV] are calculated using 
3.2 Trigger ISA JET 7.69 [28]. 

It was suggested in [20] that the triggering efficiency for events in GDM scenarios will be 
high. A full simulation of these benchmark points, which will include triggers, is under 
way [37], and we have every reason to believe that due to the generically large jet and 
lepton activity in the super symmetric events we will have trigger efficiencies in the high 
90% range for all three benchmarks. As an estimate of the trigger, in our fast simulation 
we have required each event to have at least one of the following: 
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• A jet with p T > 290 GeV, 

• three jets with px > 130 GeV, 

• four jets with px > 90 GeV, 

• one muon or stau with px > 20 GeV, or 

• two muons or staus with px > 6 GeV. 

These are the jet and muon trigger thresholds for the ATLAS detector, taken from [38], 
where we conservatively use the high-luminosity thresholds and ignore electromagnetic, 
missing-energy and tau triggers. Since one may fail to trigger on a slow stau in the normal 
running of the LHC detectors because it arrives at the trigger station too late, we addition- 
ally require triggering staus to have a velocity of /J7 > 0.9. This number is a conservative 
estimate from the ATLAS geometry, to give a trigger in the correct bunch crossing. Staus 
with ~ 1 have the longest distance to travel between RPC layers. Given a 25 ns gap 
between bunch crossings and a distance of ~ 15.5 m to the outer RPC layer the critical 
velocity is f3 = 0.67, or ^7 = 0.9, ignoring any significant energy loss. 

3.3 Background 

In order to separate the staus from background muons in the supersymmetric events, and 
to remove the Standard Model background, we use the following cuts, which we will refer 
to as the standard GDM cuts: 

• There should be exactly two stau candidates in each event, 

• these should be isolated and within the geometrical acceptance of the muon chambers 
as per the requirements for the muons listed above, 

• they should have px > 50 GeV, 

• they should have a velocity ^7 < 6.0, and 

• they should have a mass estimate from their momentum and velocity that is consistent 



with the stau mass after it has been measured (see Section 4.1 



• in addition, the sum of the number of leptons (muons or electrons) with px > 10 GeV 
and tau-tagged jets 1 with px > 20 GeV in each event should be exactly two. 

The final cut is made on the basis that in GDM scenarios with a fx NLSP, the vast majority 
of events will contain two taus in addition to the pair of staus, and these should be identified 
to fully reconstruct the particles in the supersymmetric decay chains (see Section |4~2] ) . Note 
that we have a much looser cut on /?7 than used in [36]. We still expect to remove most 
SM background because of the requirement of two stau candidates in each event, and in 
particular because of the extra requirement of two tau candidates. 



1 Assuming a 50% tau tagging efficiency. 
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After our trigger requirements and these cuts, we should have a very powerful rejection 
on Standard Model background events. We have tested this hypothesis in a fast simulation 
of a background sample consisting of the equivalent of 30 fb _1 of tt events with an NLO cross 
section of 737 pb [39], and samples of p T binned QCD, V^+jets, Z+jets and WW/WZ/ZZ 
production events. The total numbers of events are 1.75 • 10 6 for QCD and 2.5 • 10 5 for 
the other processes. We find no events that survive the cuts listed above, and, in order 
to find the first events that pass the cuts, we would need to relax the velocity cut to 
/?7 < 8.3, after which we find one WW /W Z / ZZ event that passes. However, with the 
limited statistics available, in particular for QCD, this null result must be interpreted with 
care. For example, the parametrisation of the j3 resolution assumes a Gaussian shape. If 
the true distribution has significantly larger tails, this will affect the power of the cut on f3j. 
A full simulation study of the muon velocity resolution would be required to investigate 
this possibility. It is partially because of this issue that we prefer to use a loose cut on 
/?7- Despite these problems, with the additional cuts made in the searches for various 
supersymmetric particles presented in Section |3], we expect that the cuts presented here 
will render negligible the Standard Model background in our analysis. 



4. Reconstruction of Masses in GDM Models 

In this Section, we incorporate the parametrisation of the response of the ATLAS detector 
to the slow-moving fis from Section |in a fast simulation of sparticle pair production and 
cascade decays at an LHC detector. This enables us to assess the possibility of identifying 
various heavier sparticles and the precision obtainable in mass measurements in GDM 
models. Since we have no escaping LSP in these GDM models, we are not limited to 
the standard measurements of the endpoints and shapes of kinematic distributions, see 
e.g. [40-47], and we can fully reconstruct sparticles, starting from the stau NLSP. 



We begin in Section 4.1 by looking at the starting point of all our mass measurements, 



the measurement of the stau mass. Then, in Section |4.2| we demonstrate the complete 
reconstruction of the tau momenta for events with only two taus and no hard neutrinos. 
Section 4J3 deals with the reconstruction of the decay chain q~R — > qx® — > qrfi, and 



Section 44 deals with the reconstruction of the cascade q~L — ► QX2 ~~ * 1^ ~^ Q^Xi ~~ * qUTT\. 
Finally, in Section [4.5| we discuss the possibility of measuring sneutrino and chargino masses 
using dijet decays of W bosons produced in supersymmetric cascade decays. 

4.1 Measuring the Stau Mass 

We first present an estimate of the achievable precision on the measurement of the stau 



mass for the parametrisation of the momentum and velocity resolution given in Section 3.1 



After the trigger simulation presented in Section 3.2, and making the standard GDM cuts 



of Section |3.3| with the exception of the cut on the stau mass and the requirement of two 
tau candidates, we plot in Fig. || (left) a scatter plot of the measured velocity Pj mea s versus 
the measured mass for single staus, as inferred from 

m fl = 1=-. (4.1) 



-9- 



We see a large spread in the masses for high /?7, corresponding to high momenta. This 
is expected as the assumed momentum resolution deteriorates significantly above a few 
hundred GeV, and the velocity resolution worsens at higher velocity due to the shorter ToF. 
One effect of the momentum dependence of the resolution is that the spread has a clear 
bias towards higher masses, which would be a systematic effect for a mass measurement 
using all events passing the cuts. We also see a grouping of events at high f3j and low 
masses. These are mis-measured muons that have been assigned velocities that are too 
low, and thus passed the GDM cuts used. This is particularly clear for the SM background 
events which are shown in red. 

Because the low-velocity events have a much smaller spread and less bias towards 
higher masses, it would be advantageous to use only these in a mass measurement. As we 
commented in Section 3.2, the staus from these events will be missed by the trigger system 
because they are too slow. This does not mean that they would be missed completely 
by the LHC detectors, since their muon systems are designed sufficiently robustly as to 
record hits in a large time window after a collision. However, it does mean that we will 
have to rely on other triggers, e.g., calorimetric triggers, or triggers on muons from the 
decay of the heavier sparticles. Exactly how far down in velocity it is possible to go must 
be the subject of full detector studies that lie beyond the scope of this paper. In [35] the 
velocities of the staus were found to be measurable at least down to /?7 = 0.44, albeit 
for a simulation of single staus in the muon system only. However, in [36] a conservative 
approach to the trigger issue was followed, ignoring events where the staus did not reach 
the muon chambers inside the time window of the muon triggers, and using a lower cut of 
/?7 > 0.75. Assuming that the trigger rates from our simplified trigger cuts are not entirely 
unrealistic and that we can reconstruct staus with velocities down to /?7 ~ 0.3, we make 




Figure 3: Scatter plot of measured velocity (3-f mC as versus measured mass (left), with supersym- 
mctric events in black and SM background events in red, and a corresponding plot of the measured 
stau mass (right) with an additional cut on the velocity of 0.3 < (3"f < 0.6. 
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an additional cut on (3j < 0.6 for the sample used to measure the stau mass, and plot 
the corresponding mass distribution in Fig. || (right). Fitting the mass distribution with a 
Gaussian we get a stau mass for the e benchmark of 



152.485 ± 0.021 GeV, 



(4.2) 



to be compared with the nominal value of = 152.475 GeV. 

To investigate the dependence of the measured stau mass on the lower limit on the 
velocity of staus that we can identify, we show in Fig. || the reconstructed stau mass as a 
function of the lower /?7 cut value for three different values of the upper cut on ffa(. The 
general tendency is, as could be expected, an increase in the systematic error, represented 
by the increase in distance from the nominal mass value, as the lower cut value is raised and 
low- velocity events are discarded 2 . As the lower cut approaches the upper cut value, this 
effect flattens out, while of course the statistical error increases with the reduction in the 
number of events. It is clear that, for a lower cut on /?7, in choosing the upper cut one must 
trade the loss of statistics against the reduction of systematic error from the higher- velocity 
staus. We see that there is very little change in the stau mass determination if the lower 
limit on fi'j is increased to 0.44 [35]. We also note that the error in rnf remains below the 
per-mille level even if we restrict the analysis to 0.6 < /?7 < 0.9. In all the examples shown 
in Fig. ||, the error in calculating is negligible compared with the errors encountered in 
the reconstruction of higher-mass states decaying into staus. 

Based on the clear separation of muons 
and staus in the scatter plot of Fig. [|, 
we require in the further analyses that 
stau candidates have inferred masses above 
100 GeV. This value was also used in 



r 153r 



152.9 



the background analysis of Section 3.3 



152.8- 



152.7 



152.6- 



152.5- 



152.4- 



152.3 U 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 



We use the stau mass of (4.2) in the fol- 
lowing discussions, although the much 
larger errors involved in the reconstruc- 
tion of heavier sparticles means that we 
are insensitive to the exact value, within 
the errors discussed in this Section. The 
results for the other two benchmarks 
are analogous, and the resulting masses 
can be found in Table ||. The statisti- 
cal errors are at the sub-permille level, 
which is below the expected uncertainty 
on the lepton energy scale of order 0.1%, 
and on the same level as the results 

given for various GMSB scenarios in [36]. From Fig. |], we conclude that it should be 

possible also to keep the systematic error from the momentum and velocity measurements 

at the same level if we are able to reconstruct staus with velocities down to /?7 = 0.6. 

2 It is possible that this systematic effect could be modelled and (at least partially) mitigated in a more 
complete analysis. 



Figure 4: The stau mass from a Gaussian fit as a 
function of a lower cut on (3j. The solid circles show 
the results for an upper cut of f3-f < 0.6, the triangles 
/?7 < 0.9 and the squares f3j < 1.2. The red line 
shows the nominal mass value. 
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The possibility of making further measurements on the staus by looking for staus 
trapped in the LHC detectors or in surrounding matter was discussed in [20]. Trapped, 
decaying staus could make a measurement of the stau lifetime possible, leading to an 
indirect determination of the gravitino mass, assuming the macroscopically determined 
value of the Planck scale, and could potentially even lead to a microscopic determination of 
the Planck scale and a test of supergravity if the gravitino mass could be measured directly 
from the decay kinematics [45,48]. However, these exciting possibilities lie somewhat 
outside of the scope of this paper. 



Model 


€ 


C 


V 


Mass 








h 


152.485 ±0.021 


338.24 ± 0.09 


319.02 ±0.09 




291.8 ± 1.2 


666.4 ± 1.9 


659.6 ± 1.1 


h 


307.9 ± 2.0 


674.6 ± 3.4 


673.5 ±3.5 


x\ 


176.6 ±0.9 


418.3 ± 1.8 


419.6 ±3.7 


x° 2 


339.2 ±2.0 


(796.5 ± 6.5) 


(801.0 ±7.2) 




923.4 ±3.7 






QR 


895.8 ±2.8 


1811.2 ±332.3 




b 


888.9 ±57.7 







Table 3: The mass measurements and statistical errors [in GeV] obtained in the simulations of the 
GDM benchmark points. The numbers in brackets ( ) are obtained from a fit to the edge in the 
dilepton spectrum. Jet and lepton energy scale errors are not included. 



4.2 Recalibrating Taus 



As was already mentioned in Sec- 
tion [T^, in all R-parity conserving GDM 
scenarios with a f\ NLSP, with the ex- 
ception of Drell-Yan pair productions 
of staus and some events yielding tau 
neutrinos, each event will contain two 
rs in addition to the pair of fis. In 




Section hi we discussed the measure- 
ment of the f\ mass. We can detect 
and measure the masses of many of the 
heavier sparticles by reconstructing de- 
cay chains that lead to the fis. For 
example, at benchmark point e, essen- 
tially 100% of the lightest neutralinos 
Xx decay into f\ + r, and this branching ratio is also very large at the other two points 
(~ 70%). One can reconstruct a Xx m ass peak by combining f\-r pairs, but to do this 
one needs to know the momenta of the accompanying rs, which in turn decay and lose 
momentum to escaping neutrinos. 



Figure 5: Schematic drawing of the r decay products 
Pj 1 and pj 2 , the missing momentum $t and the recon- 
structed neutrino momenta andp^, projected on 
the transverse plane. 
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Assuming that a r is relativistic in the laboratory frame, as will indeed be the case 
for any r that would pass our acceptance criteria, its hadronic decay will give a jet and a 
neutrino travelling in essentially the same direction. If there are exactly two ts in an event, 
and no other source of missing momentum, the sum of the missing transverse momenta 
from the neutrinos can be projected onto the two axes formed by the hadronic tau-decay 
jets in the transverse plane, as illustrated in Fig. [5|. Knowing the azimuthal angles of the 
two ts, one can determine the two components of each of the tau neutrino momenta in 
the transverse plane, and their momenta in the z direction follow from requiring that the 
neutrinos travel in the same direction as the jet. The same recalibration procedure can be 
used for r pair decays into one or two leptons, with leptons taking the roles of the hadronic 
jets, and the sum of two r decay neutrino momenta taking the role of the single r decay 
neutrino of the hadronic case. 

In our tau recalibration we assume that there are no other large contributions to 
the missing momentum. Possible sources is the mis-measurement of jet energies or hard 
neutrinos from the decay of heavy particles. These will change the direction and magnitude 
of the missing momentum vector in the transverse plane. This can partially be checked for 
by testing that the direction of the missing momentum lies inside the opening angle of the 
two taus. 

The result of this momentum re- 
calibration is shown in Fig. ^, where we 
plot the relative error on the tau mo- 
mentum before (dotted line) and after 
(solid line) recalibration, for the three 
benchmark points. We have used a 50% 
tau-jet tagging efficiency, and require 
that for each event the sum of the num- 
ber of tagged tau jets and leptons (no 
staus) with pj> > 10 GeV is exactly two, 
to have an unambiguous pair of tau can- 
didates for the recalibration. No trigger 
or background cuts have been used. 

We observe a dramatic improvement 
in the momentum resolution compared 
to uncalibrated taus for all three bench- 
mark points, although with a significant 
tail of taus with momenta that are too 
high. The smearing of the tau momen- 
tum is partially due to other contributions to the momentum imbalance as already men- 
tioned above, leading to a mis-calibration, but also due to the mis-tagging of tau jets and/or 
mis-identification of leptons as coming from tau decays. These mis-identifications can also 
be seen in the pre-recalibration distribution, as a tail to momenta higher than the parton 
momentum. We see a marked difference between the e benchmark and the other two, the 
e benchmark having a worse resolution. This could be expected from the relatively lower 



0.16 




(p mMS -p part )/p part 



Figure 6: The relative error between the parton- 
level and measured momenta of rs at the e (red), C 
(green) and rj (blue) benchmark points. We show the 
momenta both before (dotted line) and after (solid 
line) recalibration. 
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amount of activity in ( and rj events, where the Drell-Yan and associated production is 
more pervasive, resulting in fewer mis-identifications. 



4.3 Reconstructing the Decay q — ► qx\ — » qrfi 

The great majority of right-handed squark, qr, decays (92 %) and a few of the left-handed 
squark, <jl, decays (4 %) at the e benchmark point lead directly to the final state qx\- To 
reconstruct the masses of the lightest neutralino Xi an d the first two, almost degenerate, 
generations of squarks, we seek to isolate events with two squarks, where both decay 
according to this clean decay chain 3 

q -> qxl -> orh- (4.3) 

We first select events that survive our trigger requirements and the standard GDM cuts of 



Section |3.3| . These remove the Standard Model background, and leave us with an identified 
pair of staus. The energies of these staus are then re-calibrated to reflect the mass measured 
in Section |4.1| . Using the parametrised tau-tagging described in [34] , with an assumed tau- 
tagging efficiency of 50%, we select events where the sum of the number of tau jets with 
Pt > 20 GeV and leptons with px > 10 GeV is exactly two. These are then our basis for 



reconstructing two taus with the technique discussed in Section |4.2| . The staus and taus 
are matched according to charge, keeping only those events where there is an unambiguous 
assignment. This implies that all events with same-sign staus and all events with two tau 
jets are rejected. We do not consider any measurement of tau charge from the hadronic 
jet, leaving this as a possible improvement. However, one should note that, while this 
will certainly increase the available statistics, the lower tagging efficiency for hadronic tau 
decays compared to the efficiency of identifying leptons means that one would have roughly 
the same number of leptonic and hadronic tau candidates. Moreover, the momenta of the 
hadronic jets would be known with worse accuracy than those of the leptons. 

Calculating the invariant masses of the f-r combinations surviving these cuts we arrive 
at the distribution shown in the left plot of Fig. |7] for the e benchmark. In blue we show 
the distribution due to events with two Xi ~^ ^i + r decays, while the red distribution is 
the supersymmetric background. The total distribution shows a clear peak corresponding 
to the Xi mass. We fit the distribution by a third-degree polynomial assumption for the 
background and a Breit-Wigner distribution for the peak, giving a x? mass and statistical 
error of 

m^o = 176.6 ± 0.9 GeV, (4.4) 



which is to be compared with the nominal value of m^o = 175.2 GeV. The wide shape of 
the signal distribution is due to smearing from the momentum resolution of the staus, and 
contributes to the overestimate of the X\ mass. For the £ and r\ points the results are very 
similar, but with significantly lower statistics, even at 300 fb _1 . The masses and errors 
found are listed in Table H. 



3 We use 'clean' here in the sense that the events contain no hard leptons other than the two rs from 
the Xi decay, and only two hard jets, modulo initial-state radiation (ISR) and final-state radiation (FSR). 
They should therefore provide good measurements of the transverse missing momenta from the neutrinos 
in r decays. 
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To find the (right-handed) squark mass, we want to isolate events that contain two of 
the decays of Eq. (|4.3| ), This is done by adding to the cuts used above the following: 

• We require two jets with pt > 300, 150 GeV, 

• no other jets with px > 150 GeV, 

• missing transverse energy Et < 50 GeV, and 

• to increase statistics we relax the unambiguous assignment requirement for the r-f 
pairs, keeping those not directly excluded by charge considerations and which have 
an invariant mass within two times the peak width of the x?i F^o = 25.9 GeV. 

Here the cut on missing energy attempts to remove events with additional hard neutrinos, 
not coming from the decay of the final taus in the decay chain. 

The two hardest jets in the event are then combined with the f-r pairs and the invariant 
masses of the triples are found. The invariant mass distribution is shown on the right in 
Fig. g. The surviving events are very predominantly decays of the <jr, and the total 
distribution again shows a clear peak, now corresponding to the q~R mass 4 . Fitting the 
distribution with a second-degree polynomial for the background plus a Breit-Wigner, at 
the e benchmark point we find 

mq R = 895.8 ± 2.8 GeV. (4.5) 

This can be compared to the nominal value of rriq R = 900.8 GeV at the e benchmark 
point. The background tail to wards lower invariant masses tends to result in a small 

4 As we show below, the <jl mass can be measured using a different decay chain. 




Figure 7: The f-r (left) and jet-r-r (right) invariant mass distributions for the e benchmark point. 
For the left plot the blue distribution is for events containing two decays x? — * tti, for the right 
plot events containing two q — > qxl — > grfi decay chains. The red distributions include all other 
events surviving the cuts described in the text. 
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underestimate of the mass. We note that the accuracy of the squark mass measurement 
is dependent on the jet recalibration routine used. The statistical error after three years 
at low luminosity is lower than the systematical error expected from the uncertainty in 
the absolute jet energy scale, which is estimated to be of the order of 1% (see Chapter 12 
of [38]). For the £ benchmark point we again have a clear excess of events over background, 
and we can reconstruct a mass peak, but only containing a few dozen events. The statistical 
error for £ is large, and varies somewhat with binning and fit range. For r] only two signal 
events survive all cuts, compared to a background of ten, thus we are not able to reconstruct 
the qn in this case. We list the numerical results in Table ||[ 

By using jet-f-r combinations where the jet is tagged as coming from a b quark, we can 
estimate the b mass. Because of the low statistics obtainable with 30 fb _1 of luminosity 
we can provide only a rough estimate of a combined &i and 62 mass. Assuming a 53% 
b-tagging efficiency, with light jet rejection as discussed in Section || we have 

m h = 888.9 ± 57.7 GeV (4.6) 

We note that, whilst the production cross sections of the 61 and 62 are of the same magni- 
tude, the branching ratio of 62 - ► &X? i s m uch larger than that of bi, as 62 is mostly b^. Con- 
sequently, events with 62 dominate this mass measurement. We recall that = 883 GeV 
at the benchmark point e. For the other two benchmark points we were unable to measure 
the 6-squark masses due to the low number of events. 

Almost a half of the sparticle pair-production events at the LHC at benchmark point 
e would include gluinos, and the great majority of their decays would be into squark-quark 
pairs. One may therefore hope to reconstruct a g mass peak by plotting the invariant 
masses of jet-jet-f-r combinations. However, if we take the events lying within twice the 
width of the squark mass peak found here, and add a second jet to each jet-f-r combination 
we find that while there is a clear peak at the gluino mass in events which contain the decay 
of Eq. ( |4.3| ), initiated by a gluino, the supersymmetric background has a similar shape at 
roughly the same position and is of the same magnitude, thus preventing any gluino mass 
determination. This can be understood in terms of the softer nature of the jets from the 
gluino decay as compared to the squark decay; the background jets have a very similar 
kinematic distribution. One possible idea, which is beyond the scope of this paper to 
explore further, is the decay g — ► tt± — ► ttxi- This could be used to find the gluino mass 
if the top decays involved could be fully reconstructed. The decay has the advantage of a 
large branching ratio, but reconstructing the W involved from hadronic decays will present 
a formidable challenge. For the other two benchmark points, although the jets from gluino 
decays are somewhat harder, the background problem still remains. 

It might be possible to improve on the above results by refining the event selection, 
e.g., by using only leptonic r decays for high statistics scenarios, and/or by making a 
more sophisticated analysis of the final states, for example by modelling the expected 
supersymmetric background in a more sophisticated way than our polynomial assumptions. 
However, the results in Fig. ^ and Eqs. (4.4) and ( |4.5|) already improve significantly on those 



shown in [20], where no attempt was made to correct for the missing neutrino momenta 
in r decays, and the full kinematics of the decay was not exploited. We also observe that 
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for the e benchmark point our statistical errors are in general below the expected levels of 
errors from systematical uncertainties in lepton and jet energy scales at the LHC. 



4.4 Reconstructing the Decay q 

Although the long decay chain 



QX2 



qUXi -> qttTT\ 



qxl^qU^q^Xi^qUrn, 



(4.7) 



predominantly of the left-handed squark via a left-handed slepton, has a much lower branch- 
ing ratio than the decay of the right-handed squark discussed in the previous Section, it has 
a distinctive signature of between two and four leptons, depending on the r decay modes, 
in addition to the pair of fs. 

It is well known that a very use- 
ful step in the reconstruction of anal- 
ogous decay chains in neutralino dark 
matter scenarios is the measurement of 
the dilepton spectrum resulting from \\ - 
ttx\ decays. Accordingly, we start here 
by also looking at the dilepton edge for 
the e benchmark point, plotting the in- 
variant mass of pairs of same-flavour, 
opposite-sign leptons in the MC data. 
After the standard GDM cuts from Sec- 
tion IH to isolate events with two staus, 
loosening the cut on tau candidates to 
allow for more leptons, we pick events 
with a minimum of two additional lep- 
tons, each having px > 30 GeV. 5 Com- 
bining all pairs of these leptons with the 
same flavour and opposite sign we arrive 
at the blue distribution in Fig. ||. We 
can remove the background of uncorre- 
cted leptons by making the standard subtraction of the same flavour distribution (red) 
under the lepton universality assumption. The resulting, well known, triangular-shaped 
distribution from the decay yQ—*&£—* UX\ appears, and we can perform a fit to the end- 
point. Using a straight line fit, smeared by a Gaussian distribution to simulate sparticle 
finite-width effects and detector smearing at the edge, we measure 
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Figure 8: The different-flavour subtracted dilep- 
ton invariant mass distribution for the e benchmark 
point. We also show the opposite-sign, same-flavour 
distribution (blue, dashed line) and the opposite- 
sign, different-flavour distribution (red, dashed line). 



rn. 



118.6 ± 1.0 GeV, 



(4. 



J We have checked that the generated SM background discussed in Section 3.3 remains zero after these 
changes. 
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where m^ ax is the position of the di-lepton edge, given in terms of the masses of the 
sparticles in the decay chain by the following formula: 

(m| 9 -m})(mj - m 2 - 9 ) 
7T~ 



(mfD 2 = '—^ (4.9) 



The nominal value for the e benchmark is m^ ax = 119.7 GeV. For £ and rj the endpoints 
are found at m^ ax = 332.2 ± 9.0 GeV and m™ ax = 339.1 ± 9.8 GeV, respectively. 



With this information in hand we start the search for the decay in (4.7). First we 
attempt a reconstruction of the slepton mass in the decay £ — > £xi &t?i, using the 
standard GDM cuts of Section 13.31 with some modifications: 



• We allow for additional leptons to the tau candidates, iterating over the assigment 
of leptons as tau candidates or as "additional" , 

• we require that we have no more than two tau-tagged jets to avoid neutrinos from 
other tau decays, and 

• the tau candidates are paired with staus, and we keep only events where there is at 
least one pair with consistent charge and invariant mass within two times the width 
of the Xi mass peak, as measured in Section [4.3| . 

Following these cuts we add the additional lepton(s) to the reconstructed Xii an d show the 
resulting invariant mass distribution of the lepton-r-f triples in the left plot of Fig. ||. After 
fitting the distribution with an assumed linear background and Breit-Wigner resonance we 
find a slepton mass of 

m i = 307.9 ± 2.0 GeV, (4.10) 

which can be compared with the nominal value of = 305.9 GeV. 

The peak in the blue distribution to the left of the resonance that can be seen in 
Fig. H is due to the combination of Xi candidates with the softer lepton coming from x® 
decays, while the peak in the red distribution is due to generic soft leptons added to the 
Xi resonance. From this we suspect that a scenario with a small mj — m^o mass difference 
will encounter the same problems as the gluino mass measurement, while the x.2 mass will 
still be easily accessible. In such scenarios there will be an advantage in first measuring 
the x\ mass. 

For the other two benchmark points we get similar distributions, but again of course 
with less events and fitted masses with larger statistical errors. For £ we get = 674.6 ± 
3.4 GeV, and for rj, m- f = 673.5 ± 3.5 GeV. We note that the value for £ is quite some 
distance from the nominal value of = 683.8 GeV. The distribution of the resonance 
peak is very broad and extends towards lower invariant masses from a maximum around 
the true maximum. 

Starting from the slepton resonance we make the following additional cuts to find the 
mass of the X2 f° r the e benchmark point: 

• We take events with lepton- f-r combinations within twice the width of the slepton 
mass peak, 
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• we require events to have at least two leptons in addition to the tau candidates, 

• the additional leptons are required to have opposite sign and same flavour, and to 
have an invariant mass mu < 120 GeV, on the basis of the measured dilepton edge. 

Taking all allowed lepton pairs in an event and plotting the invariant mass of the lepton- 
lepton-f-r combinations, we arrive at the distribution in the right-hand plot of Fig. |9|. 
Fitting the clear mass peak with a Breit-Wigner function we find a m ass of 

m^o = 339.3 ± 2.0 GeV, (4.11) 

to be compared with the nominal value m^o = 339.4 GeV. 

For the £ and r\ points, the statistics are very low after all cuts. We find no real mass 
peak for £, with only a handful of events scattered over a large mass interval, and a very 
broad resonance for 77, from which one can estimate m^o = 800 ± 50 GeV. However, from 
the formula for the dilepton edge in Eq. (|4.9| ) and the x? mass found in the previous Section 
we can calculate the x% mass and statistical error for £ and 77. The results are listed in 
Table §. 

The information from the dilepton edge could also be used to restrict further the masses 
at the e benchmark by a combined \ 2 fit, but the effect is small, although the statistical 
error on the x\ mass can be halved. To increase the accuracy of the mass determination, in 
particular for the benchmark points with low cross sections, it could also be interesting to 
pursue the other invariant mass distributions that can be constructed from the leptons in 
this decay chain, e.g., the lepton-stau invariant masses and the lepton- lepton-stau invariant 
mass. A generalisation of the formulae of [47], which describe the shapes of invariant mass 




Figure 9: The lepton-f-r (left) and lepton-lepton-f-r (right) invariant-mass distributions for the 
e benchmark point. The blue distributions are for signal events, events containing a decay t — ► 
£\i — > irfi (left) and xQ t£ —* ££Xi ~ ¥ Utt\ (right), while the red distributions include all other 
events that pass the cuts given in the text. 
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distributions in cascade decays, to include massive stable end-products like the stau is also 
possible. This is, however, outside the scope of this initial study. 



What remains from the decay chain in (4.7) is to determine the left-handed squark 
mass. Taking the events in the x% mass peak, we add the two hardest jets of the event to 
the reconstructed x% Because there is a significant fraction of background events, mostly 
events where the X2 does not come from the decay of a left-handed squark, we additionally 
require a missing transverse energy Et < 50 GeV. While statistics are very low as a result of 
this selection, we find a mass peak and can perform a fit with a Breit-Wigner distribution. 



The result is plotted in Fig |10|, and we find a left-handed squark mass of 

rriq L = 923.4 ± 3.7 GeV, (4.12) 

Compared to the common mass of the first two generations of both families, mj^ ~ = 
938.2 GeV and rriu L) c L = 934.5 GeV, this fit result is slightly low, indicating that there 
may be some systematical effect, and with a relative error of the order of 1%. However, 
comparing with the squark mass found in the previous Section, this analysis clearly demon- 
strate that there are two different masses and the decay chains allow us to assign the left- 
and right-handedness as we have done. Since we could not reconstruct a X2 mass peak, 
the (Jl mass is of course also inaccessible with this technique at the £ and r/ benchmarks. 

4.5 Reconstructing v T and xf Decays 

The high fraction of the total cross section that comes from the associated production 
°f Xr~X2 P a i rs at the £ and 77 benchmark points, together with the large production of 
neutralinos and charginos in the decays of squarks at all three points, makes the dominant 
decays of X\ potentially as interesting as the x® decay via a left-handed slepton that was 




Figure 10: The jet-lcpton-lepton- t-t (left) and the W-f (right) invariant-mass distributions for 
the e benchmark point. The blue distribution is for signal events, events containing the decays of 
Eq. 4.7 and v T — > Wf\ : respectively, while the red distributions include all other events that pass 
the cuts given in the text. 
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discussed in the previous Section. Accordingly, we look now at the corresponding decays 
of the chargino. 

In the case of the e benchmark point, the decay chain X\ ~> T ^r tWt\ has a total 
branching ratio of 9.5%, and the decay x° ~~ > v T i) T — > v T Wf\, a branching ratio of 9.0%. 
This makes the search for a v T resonance feasible, if we can reconstruct hadronic decays of 
the W in the supersymmetric events. 

The reconstruction of Ws produced in the decays of heavy particles at the LHC is diffi- 
cult because of the large jet background in the high multiplicity environment. Additionally, 
since Ws will be so copiously produced in SM events at the LHC, it will be difficult to 
find effective background cuts. When a W is produced with little boost in the laboratory 
frame, the jets from the decay are relatively soft and easily drown among the other jets 
of the event. For boosted Ws the jets are harder, but have a much smaller opening angle 
in the lab frame, and are thus easily mistaken for a single jet by a cone jet algorithm. 
The quest for possible improvements in the search for hadronic decays of heavy bosons in 
supersymmetric events by the use of different jet algorithms is an interesting topic, but lies 
beyond the scope of this paper. As we shall see, because of the excellent rejection of the 
Standard Model backgrounds afforded by the properties of the stau, the identification of 
W candidates is easier in the GDM scenarios. 

To isolate events with the decay v T — > Wf\, we apply the standard GDM cuts of 
Section |3], with the exception of the requirement of two tau candidates. Instead we tighten 
the cut on stau velocity to /?7 < 2.7, which we again find removes all SM background. 
This allows us a very simple identification of events with Ws. We require one and only one 
combination of two jets in the event which has: 

• Invariant mass within two times the width, Tw = 2.12 GeV, of the W mass, m\y = 
80.42 GeV, and 

• an opening angle 9 < \ in the laboratory frame. 

The cut on angle is effective in rejecting the large background of events without a W, 
where by chance two relatively soft jets have the correct invariant mass, while, for reasons 
discussed above, it is a typical property of Ws that come from the decays of heavier 
sparticles. After these cuts and after combining the reconstructed W candidates with 



the closest stau we arrive at the invariant mass distribution shown in Fig 10. The signal 
distribution (blue) shows a clear peak above the background (red), which consists mostly 
of events without a W. By fitting the total distribution with the sum of a third-degree 
polynomial for the background and a Breit-Wigner distribution for the peak, we arrive at 
a value of the tau-sneutrino mass of 

m Pr = 291.8 ± 1.2 GeV (4.13) 

to be compared with the nominal value of m^ T = 288.5 GeV. The overestimated value for 
e is due to the long tails of the signal distribution to the right of the peak, which in turn 
comes from combinatorial effects, such as pairing the W with the wrong stau candidate. 
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From events in the tau-sneutrino resonance peak we should in principle be able to 
reconstruct the chargino mass by adding tau jets to the reconstructed sneutrinos. However, 
the large background and the soft nature of these taus means that the signal drowns in a 
background of the same shape and position, much as in the case of the gluinos. For the £ 
and r\ benchmark points the situation is similar, we can reconstruct a v T mass peak, but 
no chargino. The sneutrino masses found for these points are listed in Table |3[ 



5. Conclusions 

We have demonstrated in this paper that the LHC detectors (as exemplified here by the 
ATLAS detector), would have the abilities to make many valuable measurements in bench- 
mark supersymmetric scenarios with gravitino dark matter and a f NLSP. The triggers 
planned would record GDM events with high redundancy and efficiency, and the robust 
designs of the detectors would enable them to characterise very well fs produced in LHC 
collisions. Once recorded, the f tracks in GDM events could be identified with high ef- 
ficiency, and useful momentum, ToF and dE/dx measurements would be made. We find 
that it should be possible to measure the stau mass with statistical and systematic errors 
at the sub-permille level, assuming only the reconstruction of staus with velocity down to 
/?7 rj 0.6. The majority of GDM events with a pair of f NLSPs would also contain a pair 
of r leptons, which should be detectable with good efficiency at the LHC. Combining these 
with the f measurements would enable first the Xi an d subsequently heavier sparticles 
to be reconstructed and their masses measured with accuracies that often approach the 
expected systematic energy scale uncertainty of jets and leptons. We demonstrate these 
abilities with the aid of r momentum recalibration using transverse momentum balance, 
and making appropriate kinematic cuts on the final-state kinematic distributions in spar- 
ticle cascade decays. Our results are consistent with the preliminary estimates of sparticle 
observability made in [20]. 

The LHC will be the first accelerator to explore a new energy domain, and there are 
high hopes that it will reveal new physics beyond the Standard Model. Supersymmetry 
is among the most 'expected' new physics surprises, which has been much studied by 
the LHC collaborations as well as by theorists. Most of these studies have assumed the 
'expected' missing-energy signature of R-conserving supersymmetry with a neutralino LSP, 
though other possibilities have also been explored. It has been realised recently that some 
of these other possibilities may be as generic as missing energy, and the scenario studied 
here of gravitino dark matter and a stau NLSP is just one of many possible 'unexpected' 
possibilities. Nevertheless, it is encouraging that, although the LHC detectors were not 
designed with such an exotic scenario explicitly in mind, they seem to be capable also 
of dealing with such an 'unexpected' physics surprise. We have shown here that valuable 
insight into such a gravitino dark matter scenario could be obtained already with a relatively 
modest initial amount of LHC luminosity, and so could inform the design of possible future 
experiments at the high-energy frontier. 
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